Thermodynamics of spin-orbit-coupled Bose-Einstein condensates 
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In this paper we develop a quantum field approach to reveal the thermodynamic properties of the 
trapped BEC with the equal Rashba and Dresselhaus spin-orbit couplings. In the experimentally- 
feasible regime, the phase transition from the separate phase to the single minimum phase can be 
well driven by the tunable temperature. Moreover, the critical temperature, which is independent 
of the trapped potential, can be derived exactly. At the critical point, the specific heat has a large 
jump and can be thus regarded as a promising candidate to detect this temperature-driven phase 
transition. In addition, we obtain the analytical expressions for the specific heat and the entropy 
in the different phases. In the single minimum phase, the specific heat as well as the entropy are 
governed only by the Rabi frequency. However, in the separate phase with lower temperature, we 
find that they are determined only by the strength of spin-orbit coupling. Finally, the effect of the 
effective atom interaction is also addressed. In the separate phase, this effective atom interaction 
affects dramatically on the critical temperature and the corresponding thermodynamic properties. 

PACS numbers: 03.75.Mn, 03.75.Hh, 67.85.-d 



I. INTRODUCTION 



The spin orbit coupling (SOC), which describes the in- 
teraction between the spin and orbit degrees of freedom of 
a particle, has not only generated many interesting quan- 
tum phenomena in modern physics ranging from the nu- 
clear physics to condensed-matter physics, and but also 
become an important resource for realizing fault-tolerant 
topological quantum computing By controlling the 
external lasers, the different kinds of SOCs have been 
proposed to be simulated in the trapped Bose-Einstein 
condensates (BECs) with the neutral atoms @. Espe- 
cially, in recent experiment at NIST, the equal Rashba 
and Dresselhaus SOCs has been realized successfully in 
the ultracold 87 Rb atoms by a couple of Raman lasers 
0. Attributed to this pioneer experiment, the investi- 
gation of SOC-driven BECs has attracted much atten- 
tions. Moreover, rich many-body phenomena with no 
analogy in condensed-matter physics (in BECs, all ultra- 
cold atoms can occupy the same quantum state) have 
been predicted by considering the ground-state proper- 
ties [H-I23I]. For example, in the presence of the equal 
Rashba and Dresselhaus SOCs, the BEC is made up of 
two non-orthogonal dressed atom spin states carrying dif- 
ferent momenta. Furthermore, the interaction between 
these spin states are modified, driving a quantum phase 
transition from a spin-mixed state to a phase-separated 
state 24|. In fact, even if the effective atom interaction 
governed by both the inter- and intra- spin interactions 
is not taken into account, a quantum phase transition 
from a separate phase (SP) to a single minimum phase 
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(SMP) can also occur [25|,|26|. In very recent experiment, 
this new quantum phase transition has been observed by 
measuring the amplitude ratio of spin and momentum 
oscillation [27j . 

It has been known that quantum phase transitions gov- 
erned by the ground-state energies occur at absolute zero 
temperature [28[. However, it is unattainable experimen- 
tally due to the third law of thermodynamics, i.e., any 
system must work at a finite temperature. Thus, it is cru- 
cially important to investigate the thermodynamic prop- 
erties to fully understand the fundamental physics for a 
given system. For instance, in the framework of finite- 
temperature theory, the system's real evaluation can be 
described more accurately and some important physical 
quantities such as the specific heat, the entropy and the 
free energy, which have no zero temperature correspon- 
dence, can be explored. More importantly, some exotic 
phenomena driven only by thermal fluctuations can be 
revealed 

Motivated by the experimental developments and the 
third law of thermodynamics, we, for the first time, 
develop a quantum field approach to reveal the ther- 
modynamic properties of the trapped BEC with the 
equal Rashba and Dresselhaus SOCs. Our main results 
are given as follows: (I) In the experimentally-feasible 
regime, the phase transition from the SP to the SMP can 
be driven by the tunable temperature. Moreover, the cor- 
responding critical temperature is derived exactly and is 
independent of the trapped potential. (II) We find that 
the specific heat has a large jump at the critical temper- 
ature. This step behavior is quite different from that of 
the atom population, which varies smoothly when cross- 
ing the critical point. It implies that the temperature- 
driven phase transition can be well detected by measur- 
ing the specific heat. (Ill) In the different phases, the 
analytical expressions for the specific heat and the en- 
tropy are also given. In the SMP, the specific heat as 
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FIG. 1: (Color online) (a) The experiment setup for realizing 
the equal Rashba and Dresselhaus SOCs in the trapped BEC 
at NIST ||. (b) The energy level structure of 87 Rb atoms. 

well as the entropy are governed only by the Rabi fre- 
quency However, in the SP, the strong SOC modifies 
the energy structure and thus the thermodynamic statis- 
tics. At lower temperature, we find that the specific heat 
and the entropy in such phase are determined only by 
the SOC strength. (IV) Finally, the effect of the effec- 
tive atom interaction is also addressed. In the SMP, no 
collective excitations can be found in SOC-driven BEC 
and thus the effective atom interaction does not affect 
the thermodynamic properties. However, in the SP with 
SOC-induced macroscopic excitations, this effective atom 
interaction affects dramatically on the critical tempera- 
ture as well as the other thermodynamic quantities. For 
example, for the repulsive atom interaction, the critical 
temperature decreases, and vice versa. 

II. MODEL AND HAMILTONIAN 

Figure 1 shows the experimental scheme about how 
to create SOC in the trapped BEC with the ultracold 
87 Rb atoms at NIST 0. In their experiment, the BEC 
is trapped in the xy plane through a strong confinement 
with frequency uj z along the z direction. In the large 
detuning A, the momentum-sensitive coupling between 
two hyperfine ground states \F — \,mp = — l)(|f)) and 
\F = \,mp = 0)(||)) is constructed by a pair of Ra- 
man lasers with Rabi frequencies fii and f2 2 incident 
at a 7r/4 angle from the x axis, as illustrated in Fig. 
1(a). In the dressed-state basis |f) = exp (iki • r) |f) and 
1 1) = exp (ik 2 • r) \\), where ki and k 2 are the wavevec- 
tors of the Raman lasers, an effective SOC, which is iden- 
tical to the one-dimensional equal Rashba and Dressel- 
haus SOCs in condensed-matter physics, can be achieved. 



Moreover, the corresponding Hamiltonian with the atom- 
atom collision interaction can be written from the cou- 
pled Gross-Pitaevskii equations as [25| 

Ho = hui x Na) a + hQS x — 70 V mftw x i{a* — a)S z H — -S z . 

(1) 

H ere, a^a is a harmonic trap mode with a = 
\J mui x /2h{x + ip x /mui x ) and m being the atom mass. 
S z = (<f>!$t - <f>|$j,)/2 reflects the experimentally- 
measurable population between the different spin com- 
ponents. uj x is the trapped frequency in the x direc- 
tion, fl = fiif^/A is the effective Rabi frequency. 70 = 
\f2hkpjra with hk^ = ^/2i:h/X being the SOC strength, 
where A is the wavelength of the Raman laser. The effec- 
tive atom interaction q is proportional to N(g^ + gn — 
2g n ), where g n = g n = 4irh 2 N(c Q + c 2 )/(ma z ) and 
.9-1-4- = 47rfi, 2 A r co/ma 2 are the inter- and intra- spin inter- 
action constants with cq and c 2 being the s-wave scatter- 
ing lengths and a z = \J1-KhjmLo z . N is the total atom 
number. 

If defining the number-dependent trapped frequency 
lu = Nlo x and the effective SOC strength 7 = y / ro7o, 
Hamiltonian ([T]) can be rewritten in the rotating frame 
as (H = 1 henceforth) 

H = uja)a + nS z + ^^{a) + a)S x + ^-St (2) 

In the following discussion, we focus mainly on Hamilto- 
nian ([2]), in which (S x ) stands for the atom population. 
Before proceeding, we estimate the relative parameters 
under current experimental conditions [1, H3, [3(| . In the 
experiment of NIST, the tunable trapped frequency u> x is 
of the order of 10 Hz, and correspondingly, the number- 
dependent trapped frequency to is of the order of MHz 
for N = 1.8 x 10 5 . Parameter 7 2 is of the order of 
kHz for A = 804.1 nm. The effective Rabi frequency 
Q can range from zero to the order of MHz. In addition, 
since Co = 100.86 as and c 2 = —0.46 as with as be- 
ing the Bohr radius, we have g^ ~ g^ ~ g-^ and thus 
q ~ 0. It means that the effective atom interaction need 
not be taken into account in the NIST's experiment. It 
should be pointed out that this effective atom interac- 
tion can be well controlled through Feshbach resonance 
[3lj . Moreover, its magnitude can reach the order of 
MHz near the Feshbach resonant point. Finally, we will 
take El = h 2 kj^/2m, which is of the order of kHz, as the 
natural unit of the energy for simplicity. 

III. THERMODYNAMIC EQUILIBRIUM 
EQUATION 

A key step to extract the thermodynamic properties of 
the SOC-driven BECs is to obtain the partition function 
of Hamiltonian ^ [29| . Here we develop a quantum 
field approach, i.e., an imaginary-time (r = it) func- 
tional path-integral technique, to arrive at the target. 
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We first rewrite the collective spin operators in the rep- 
resentation of the Grassmann Fermi fields, namely, S z = 



where the Fermi operators £tj(^;) and ^(z^i) satisfy the /-/3 w 

anticommutator relations {/ij, H } = {v\, Vj ) = 8 ly Fur- A = J dr ^* d ^ + Y,^ d ^* + v * d r v i) + h f], (4) 
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In Eq. (HJ, [d7?(r)] = tp* , fi, n* , v, v*] is the path 
integral measure. The Euclidean action is given by 



and S- = S], 



thermorc, we transform the harmonic trap mode a^(a) 
into a single mode bosonic field As a consequence, 

the partition function is obtained by 



[dr){T)} exp [-A(t)] 



(3) 



where d T = d/dr, (3 = l/(fcsT) with ks being the Boltz- 
mann constant and T being the system's temperature, 
and 



H, 



LJ1p*1p + Q } (11*11 



N 
i=l 



N j- 

2— 1 V 



(5) 



r 



Since Hamiltonian ([2]) has two degrees of freedom includ- 
ing the spin and orbit cases, it is very difficult to directly 
discuss the partition function to extract its fundamental 
thermodynamic properties. The usual method is that we 
eliminate one de gree of freedom by integrating the Eu- 
clidean action A [29|]. Without the effective atom inter- 
action (q = 0), the Euclidean action A is a quadric term 
and the corresponding integral is Gaussian. It means 



that in this case we can integrate over the Grassmann 
Fermi fields and then obtain the partition function of the 
bosonic mode. However, for nonzero q (q ^ 0), the inte- 
gral in the Euclidean action A is not Gaussian and thus 
the corresponding integral is hard to be solved directly. 
Here we introduce an auxiliary field x to circumvent this 
difficult. Based on this auxiliary field x, we have [32| 



exp[--^ Yl ~ v * v if \ K J \ dr i\ exp{ J dr[^-x*x - 



1 N 



(6) 



In analogy of the mean field approximation, the value 
of auxiliary field x determines (S x ), as will be shown. 
Substituting the formula about the auxiliary field x into 
the Euclidean action A yields 

rP 

A(j>,x) = A (iP,x)+ V / dT$*G(<ip,x)$i, (7) 
where 



with 



A (tp,x)= / dr [ip* (u + d T ) ip — x*x] 
Jo 



and 



G(^,x) 



T(ip,x) dr-il 



(8) 



(9) 



(10) 



ry>,x) = i x i-y>* + ii>)- 



(x*+x). (11) 



For the effective Euclidean action A in Eq. ([7]), we 
can integrate over the Grassmann Fermi fields, i.e., the 
degree of freedom for the spin, and then obtain 

A = N dr [** (uj + d T ) * - X*X - TrlnG] , (12) 
Jo 



where = ip/y/N and X = x/vN. Finally, by means 
of the standard stationary phase approximation, namely, 
6A/6V = SA/S^* = and SA/SX = SA/SX* = 0, the 
required and X, which play a crucial role in determin- 
ing thermodynamic properties of Hamiltonian ([2]), can 
be obtained by 



v|> = <3>* = 

X = X* = 



|( 7 2 *-^X)tanh(f) 
3(7«Vw*-gJC)tanh(^) 



(13) 
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FIG. 2: (Color online) The scaled atom population \{S X )\/N 
as the functions of the effective SOC strength 7 and the tem- 
perature T without the effective atom interaction (q = 0), 
where the effective trapped frequency uj ~ 5.1 x 10 3 .El and 
the effective Rabi frequency f2 = 5.0i?L- The red dashed line 
given by Eq. (|f 6|) determines the critical boundary. Inset: 
The scaled atom population \{S X )\/N as a function of the 
temperature T with the effective SOC strength j 2 = 2.6-El- 

where 

C = \/si 2 + 4( 7 V^f- A) 2 . (14) 

It should be noticed that in the derivation of Eq. (| 13[) we 
focus on the constant path that ^ is not influenced by r, 
namely, d T ^f = [2{|. According to Eq. (|T5|) we have 

._2&Lj> (15) 

Equation ([15)) shows clearly that there exist a trivial 
solution = *]/* = 0, and the nontrivial solutions 
\J> = \J/* = ±\£ and A = q^/u^o/j when 7 / 0. 
Moreover, these nontrivial solutions are governed by the 
nonlinear equation Co/ [2 ( 7 2 — q)] = tanh(/3£o/2), where 

Co = v /fi2 + 4 (7v / w*o-^o) 2 - With the help of the sta- 
ble condition at the equilibrium points, we can obtain 
the required solutions of both ^ and X and thus reveal 
the thermodynamics of Hamiltonian [29[ . 

IV. WITHOUT EFFECTIVE ATOM 
INTERACTION 

We first address the case of q = 0, which has been 
realized at NIST ||. At zero temperature (T = 0), 
tanh(/?C/2) = 1 and thus, Eq. (fl"5"j) becomes <!> = 
27 2 \1//C, which leads to solutions of * = (S x ) = 
for 7 < ^fV2 and * = ±^(4 7 4 - ft 2 )/(4w7 2 ) and 
(&,) = -^1 -ft 2 /(4 7 4 ) for 7 > These zero- 

temperature solutions agree well with the direct numeri- 
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FIG. 3: (Color online) The specific heat C v (Red lines) as 
well as the entropy S (Blue dashed lines) as a function of 
the temperature T without the effective atom interaction 
(q = 0). The plotted parameters are given by ui ~ 5.1 x 10 3 -El, 
fl — 3.0-El, and j 2 = 1.8-El. These lines and the circles 
(squares) denote the numerical simulation (NS) and the ana- 
lytical result (AR), respectively. Inset: The first-order deriva- 
tive of S/Nks versus the temperature T (Blue dash-dotted 
line). 



cal simulation of the SOC-driven Gross-Pitaevskii equa- 
tions [25[ . The nontrivial variations of both \1/ (atom mo- 
mentum) and (S x ) (atom population) show that a quan- 
tum phase transition occurs by adjusting the effective 
SOC strength 7. Moreover, we can call \& = (S x ) = 
as the single minimum phase (SMP) with no collective 
excitations, and "J ^ and (S x ) ^ as the separate 
phase (SP) with the macroscopic excitations (2||. With 
the increasing of the temperature T, the order parameter 
\& or (S x ) will be destroyed by thermal fluctuation. In 
particular, when ^(Tc) = (S x ) (T c ) = 0, the system en- 
ters into the SMP from the SP. By means of *(T C ) = 0, 
the critical temperature can be obtained exactly by 



2k B arctanh(2^r) ' 

Eq. (fTB|) shows that the phase transition from the SP 
to the SMP can be driven by the tunable tempera- 
ture. Moreover, the corresponding critical temperature 
obtained exactly is independent of the trapped potential 
uj. When SI = 0.2£"l and 7 2 = El, the critical temper- 
ature is evaluated as T c = 84.9 nK, which is feasible in 
experiments about SOC-driven BECs. 

Having obtaining the critical temperature, we discuss 
the experimentally-measurable atom population at finite 
temperature. In terms of the thermodynamic equilibrium 
equation (|15|) . the partition function is given in the SMP 
with #(T) = by 

Z SM p = cxp{-iV0[-|ln(2cosh(^))]}, (17) 
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FIG. 4: (Color online) The critical temperature T c as a 
function of the effective SOC strength 7 for the different 
effective atom interactions q — — 0.3i?L (Black solid line), 
q = 0.0-El (Red dashed line) and q = 0.3-El (Blue dash- 
dotted line). Inset: The critical temperature versus the ef- 
fective atom interaction q with the different effective SOC 
strength -y 2 = 2.6E L (Black solid line) and j 2 = 5.0£ L (Red 
dotted line). In these figures, the other plotted parameters 
are given by uj ~ 5.1 x W 3 E L and Q = 5.0E L . 



{S x f P (T) = -^*(T) 



(21) 



in the SP. In general, 'F(T') shoud be determined numer- 
ically by solving the nonlinear equation (|15[) . However, 
when T = 0, Eas . (T2 0T) and (f2Tj) reduce to the known 
analytical results [25|, [26|. In Fig. ® we plot the scaled 
atom population \(S X )\/N as the functions of the effec- 
tive SOC strength 7 and the temperature T. This figure 
shows that thermal fluctuations destroy the collective ex- 
citations. As a result, the system finally enters into the 
SMP from the SP. 

For a full understanding of the temperature-driven 
phase transition, it is very important to discuss the ther- 
modynamic quantities in the different phases. Here we 
consider the specific heat per atom and the entropy per 
atom. The other thermodynamic quantities can be calcu- 
lated using the same procedure. By means of the formula 



whereas it becomes 



Z SP = 2cxp{-iV/3[c^--ln(2cosh(^-))]} (18) 

in the SP with ^>(T) ^ 0. Thus, the atom population 
can be derived from the formula 



dQnZ) 



^N/3d(2 lv ^^) 



by 



(S X ) SMP (T) = 



(19) 



(20) 



with U 



(22) 



\nZ being the total energy, the specific 



heat per atom in the SMP is obtained exactly by 



/-rSMP 



o 2 



2fc s T 2 



sech 



2k B T 



(23) 



which is independent of both the trapped frequency ui 
and the effective SOC strength 7. In the SP, the specific 
heat per atom is evaluated as 



1 



2fc B T 2 



c + 



sech 2 



2k B T 



C" 
kuT 



2C ] tanh 



2k B T 



- Auj 



k R T 



~k^f 



where C = \f^ 2 + 4tJ7 2 * 2 , (' = d(/df3 = 47 2 w**'/C, as expected, 
and C" = d 2 C/d(3 2 = 4 7 2 w (il 2 ^' 2 + C 2 **") /C 3 with 



*' = dV/dp = C 2 /{2w* [l- 2/37 2 + cosh(/3C)]} and 
tf" = d 2 ^/d(3 2 = {-16 7 4 o; 2 [l-/?7 2 sech 2 (^C/2)]* 2 *' 2 + 
7 2 Csech 2 (/3C/2)tanh(^C/2)(C 2 + 4/3 7 2 w1"I'') 2 + 
4 7 2 wC 2 [l ^ 7 2 sech 2 (/3C/2)]*' 2 
87 4 wC 2 sech 2 (^C/2)**'}/{-4 7 2 wC 2 [l 
/37 2 sech 2 (/?C/2)]*}. The specific heat C v p implies 
that the strong SOC can modify the energy structure of 
Hamiltonian and thus the thermodynamic statistics, 



It is very hard to directly extract the fundamental 
properties of the specific heat Cy from the complicate 
expression (j2~4"|) . In Fig. [3J we plot the specific heat Cy 
(Red lines) as a function of the temperature T. This fig- 
ure shows two interesting features. (I) The specific heat 
Cy has a large jump at the critical point T c , separating 
the SP from the SMP. This step behavior is quite different 
from that of the atom population (S x ) in Fig. [5J which 
varies smoothly when crossing the critical point. It im- 
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FIG. 5: (Color online) The scaled atom population K-S^)! /N 
as a function of the temperature T for the different effective 
atom interactions q — — 0.3-Bl (Black solid line), q = 0.0i?L 
(Red dashed line) and q = 0.3-Bl (Blue dash-dotted line). 
The other plotted parameters are given by lo ~ 5.1 x I0 3 i?L, 
fl = 5.0-El and y 2 — 2.6-El, respectively. 



Thus, the specific heat in the SP can be obtained ana- 
lytically by 



/-rSP 



27 4 
k R T 2 



sech 



7 



k R T 



(25) 



which agrees well with the direct numerical simulations, 
as shown the Red lines of Fig. [3] Eq. (|25[) shows, in 
contrast to the behavior of the specific heat C V MP , the 
specific heat C v p is governed only by the effective SOC 
strength 7 , i.e., it is independent of both the trapped 
potential uj and the effective Rabi frequency Vt. 

Another important thermodynamic quantity discussed 
in this paper is the entropy, which can be derived from 
the formula 



S 



dG 
df 



(26) 



with G = — k B T hi Z being the Gibbs function. In the 
SMP, the entropy per atom is obtained exactly by 



plies that the temperature-driven phase transition can be 
well detected by measuring the specific heat Cy- (II) At 
lower temperature, we have approximately ~ , F" ~ 0. 



S SMP = 2k R ln 



2 cosh 



n 



2k B T 



whereas it becomes 



■ tanh 



T \2k B T J ' 
(27) 



S 



SP 



In 2— + 2k B In 

N 



2 cosh 



2k B T 



in the SP. At lower temperature, the entropy can be eval- 
uated approximately as 



?sp. 



fin 2 
i2kB <2N~ 



-In 



2 coshl 



k B T 



k B T 



-tanhl 



knT 



(29) 

Similar to the behaviors of the specific heat Cy, the en- 
tropy 5* in the SMP is governed only by the effective Rabi 
frequency f2, whereas it is determined only by the effec- 
tive SOC strength 7 in the SP with lower temperature. 
However, its step behavior at the critical point T c is very 
small, as shown the Blue lines of Fig. [31 



—tanh ( - 
T \2k B T 

I 



c+ 



c 



knT k R T 2 



(28) 



In Fig. 2] we plot the critical temperature T c as a func- 
tion of the effective SOC strength 7 for the different ef- 
fective atom interactions. This figure shows clearly that 
for the attractive interaction (q < 0), the critical tem- 
perature T c increases for a fixed SOC strength 7, and 
vice versa. It means that in experiments we can ma- 
nipulate the effective atom interaction q to arrive at the 
experimentally- required critical temperature T c , where 
the thermodynamic phase transition from the SP to the 
SMP occurs. 



V. WITH EFFECTIVE ATOM INTERACTION 

In this section, we illustrate the effect induced by the 
effective atom interaction, which is indeed controlled by 
the Feshbach resonant technique in experiments. With 
the effective atom interaction (q ^ 0), the critical tem- 
perature can be also obtained exactly by 



Tc(q) 



O 



2k B arctanh[ 2 ^ 2 _ 



(30) 



On the other hand, in the SMP, no collective excita- 
tions can be found in SOC-driven BEC and thus the effec- 
tive atom interaction does not change the energy struc- 
ture. It implies that the partition function in the SMP 
is the same as Eq. (fTT|) , i.e., it is independent of the 
effective SOC strength 7, the trapped frequency uj and 
the effective atom interaction q. However, in the SP the 
strong SOC leads to the system's collective excitations 
with nonzero atom population. As a result, the term 
jfS 2 in Hamiltonian plays a crucial role in system's 
energy structure and thus the thermodynamic statistics. 
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FIG. 6: (Color online) The specific heat CV as a function of 
the temperature T for the different effective atom interactions 
q = -0.25£l (Black solid line), q = 0.00£l (Red dashed line) 
and q = 0.25-El (Blue dash-dotted line). The other plotted 
parameters are given by co ~ 5.1 x 10 3 i?L, = S.OE'l and 
7 2 = 1.8-El, respectively. The three gray lines indicate the 
critical temperatures in the different cases. 



In such case, the partition function becomes 
Z SP (q) = 2exp{-iV/3[^ - % ln[2 cosbA]]}, (31) 



where 6(q) = j 2 — q, $ and £ can be obtained from the 
nonlinear equation (|15j) with the effective atom interac- 
tion q. 

Based on the obtained partition function, the 
experimentally-measurable atom population is given by 
(S x f MP (q,T) = and (S x f P (q,T) = -^(q,T), 
which are plotted in Fig. [SJ It is clearly that in the 
presence of the attractive interaction (q < 0), the critical 
temperature increases and the system is more inclined 
to locate at the SP, and vice versa. This conclusion 
is identical to the result of Fig. @] (along the gray 
dotted line j 2 = 2.6-El there). In addition, in the 
SMP, the specific heat Cy and the entropy S are also 
identical to the results of q = 0. Whereas, in the SP 
they become C v p (g) = ^f-[(( + PC'f scch 2 (j3(/2) + 
2 (2C + /3C") ta nh(/3C/2)-4w (l - q/-/ 2 ) (2**' + /3*' 2 + 
/3tftf")] and S SP (q) = fc s {ln2/iV + 21n[2cosh(/3C/2)] 



/?(C + /?C')tanh(/3C/2) 

C(q) 
*'(«) 



2^ 2 w(l - g/7 2 )*^'}, where 



= 7 2 C 2 /{2<5a; [1 - 2/35 + cosh (/?£)] and 



V"{q) = {-16?7 4 a; 2 [l - /35sech 2 (/3C/2)]* 



<5Csech 2 (/3C/2) tanh (/?C/2) [C 2 



4^77 2 w** 
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/3&ech 2 (0C/2)]*' 2 



4t7 2 wC 2 [1 

Sj^uiChectf (/3C/2) **'}/{-4wr ? 2 C 2 [l 
/3<5sech 2 (X/2)]*} with 77(g) = <5(g)/ 7 
When T < T c , they reduce to the forms 



r<SP 



(?) 



2<P 



sech 2 (A) . 



7 - 9/7- 



(32) 



B \ ^ + In [2 cosh (A)] - A tanh (A) 



(33) 

where A(g) = 6(q)/kBT. In Fig. HI the specific heat 
Cy as a function of the temperature T for the different 
effective atom interactions is plotted. This figure shows 
that the fundamental properties of temperature-driven 
phase transition remain in the framework of the effec- 
tive atom interaction. However, for the repulsive inter- 
action {q > 0), the critical temperature T c can decrease. 
Moreover, the step amplitude at the critical point also 
decreases. For the attractive interaction (q < 0), the op- 
posite results exist. The similar behaviors of the entropy 
S can be also found. 



VI. CONCLUSIONS AND REMARKS 

In summary, we have explored the thermodynamic 
properties of the trapped BEC with the equal Rashba 
and Drcssclhaus SOCs, which has been realized in ex- 
periments. The thermodynamic phase transition from 
the SP to the SMP as well as the critical temperature 
has been revealed. We have also discussed the important 
thermodynamic quantities such as the specific heat and 
the entropy and obtained their analytical expressions in 
the different phases. At the critical point, the specific 
heat has a large jump and can be thus regarded as a 
promising physical quantity to detect this temperature- 
driven phase transition. Finally, we have illustrated the 
effect of the effective atom interaction, which can be well 
controlled by the experimentally-feasible Feshbach reso- 
nant technique. Especially, we have found that in the 
SP this effective atom interaction affects dramatically on 
the critical temperature and the corresponding thermo- 
dynamic properties for the SOC-driven BEC. Before end- 
ing up this paper, we briefly make two remarks. Firstly, 
our analysis is mainly based on the imaginary-time func- 
tional path-integral approach, in which the fluctuation of 
the space-dependent physical quantity (such as the den- 
sity) is usually "hidden" or "averaged out" |2{|. As a 
consequence, the important stripe phase predicted be- 
fore 0, HH, cannot be distinguished from the SP ef- 
fectively. Secondly, without SOC (7 = 0), Hamiltonian 
(|2|) turns into H = HS Z + jfS 2 . In this case, the rela- 
tion Xq — q^/uj^o/j becomes invalid. However, we can 
use the same procedure (introducing the auxiliary field 
in path-integral technique) to discuss the corresponding 
thermodynamics. 
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